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Abstract. Several recent observations using standard rulers and standard candles now suggest, either individually 
or in combination, that the Universe is close to flat, i.e. that the curvature radius is about as large as the horizon 
radius (~ 10h~^ Gpc) or larger. Here, a method of distinguishing an almost flat universe from a precisely flat 
universe using a single observational data set, without using any microwave background information, is presented. 
The method (i) assumes that a standard ruler should have no preferred orientation (radial versus tangential) to 
the observer, and (ii) requires that the (comoving) length of the standard ruler be known independently (e.g. 
from low redshift estimates). The claimed feature at fixed comoving length in the power spectrum of density 
perturbations, detected among quasars, Lyman break galaxies or other high redshift objects, would provide an 
adequate standard candle to prove that the Universe is curved, if indeed it is curved. For example, a combined 
intrinsic and measurement uncertainty of 1% in the length of the standard ruler L apphed at a redshift of z = 3 
would distinguish an hyperbolic (fim = 0.2, SIa ~ 0.7) or a spherical {Qm = 0.4, Qa = 0.7) universe from a flat 
one to 1 — P > 95% confldence. 
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1. Introduction 

In the Friedmann-Lemaitre-Robertson- Walker model 
(Weinberg 1972), the Universe is an almost homogeneous 
3-manifold of constant curvature. This manifold may be 
the 3-hyperboloid H'^, flat Euclidean space or the 3- 
sphere S^, or a quotient manifold of one of these, e.g. the 
3-torus T3 (Schwarzschild 1900, 1998). 

A directly geometrical way to measure curvature is by 
a standard candle or a standard ruler, i.e. a class of ob- 
jects of which the intrinsic brightness or comoving length 
scale is believed to be fixed. Several recent applications of 
standard candles or standard rulers include the standard 
candle applications of Perlmutter et al. (1999) and Riess 
et al. (1998), the standard ruler applications of Roukema 
& Mamon (2000, 2001), Lange et al. (2000) and Balbi et 
al. (2000), which use geometrical information in the tan- 
gential direction, and the standard ruler application of 
Broadhurst & Jaffe (2000), which uses geometrical infor- 
mation in the radial direction. 

These recent measurements, individually or in com- 
bination, favour an 'almost' flat local universe. However, 
J&Se et al. (2000) reject a flat universe to just under 95% 



significance, finding Q 



tot 



1.11 



+0.13 



('95% confidence'). 



Whether the Jaffe et al. (2000) result is just due to 
random or systematic error and the observable Universe is 
in fact flat to high precision as predicted by many models 
of inflation, or whether the Universe really is measurably 
curved, distinguishing an almost (but not) flat universe 
from an 'exactly' flat model will require considerably more 
precise techniques than have been previously applied. 

Given the claimed existence of a comoving standard 
ruler (Broadhurst & Jaffe 2000; Roukema & Mamon 2000, 
2001), what could possibly be the most model free tech- 
nique for testing the flatness hypothesis is presented here, 
for the case where the Universe is curved and indepen- 
dence of orientation of a comoving standard ruler is used 
to refute the flat universe hypothesis. Since 

(i) a standard ruler should have no preferred orientation 
(radial versus tangential) with respect to the observer, 
and 

(ii) the radial and tangential comoving distances differ by 
a sin or sinh factor if the Universe is curved [eq. (5)], 

then if the lengths of a standard ruler in the tangential 
and radial directions are proved to be unequal under the 
assumption of a flat universe, it might be thought that 
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this would falsify the flat universe hypothesis. However, 
as shown below, the additional constraint that 

(iii) the standard ruler should have a known, fixed comov- 
ing value 

is required in order to falsify the flat universe hypothesis. 

Since the claimed standard ruler is in the linear regime 
of density perturbations, and is presumably a primordial 
feature in the power spectrum (if real), its inability to 
evolve in comoving length scale in less than a Hubble time 
would make it free of the evolutionary effects present for 
standard rulers or standard candles defined by collapsed 
objects. 

Alcock k Paczyiiski (1979), Phillipps (1994), 
Matsubara & Suto (1996) and Ballinger, Peacock & 
Heavens (1996) have previously pointed out the po- 
tential usefulness of (i) and (ii), and have suggested 
applying these at quasi-linear or non-linear scales, i.e. at 
r ^10h~^ Mpc, including some analysis of how to try to 
separate out peculiar velocity effects. However, they did 
not discuss how to lift the degeneracy in the two curvature 
parameters (flnnf^A) which remains after using (i) and 
(ii), and the problem of evolution of the r < 10 Mpc 
auto-correlation functions of galaxies and quasars offers 
potentially serious systematic uncertainties. 

In contrast, by using a comoving standard ruler in the 
linear regime, the constancy in the comoving scale over a 
Hubble time [(iii) above] provides an additional constraint 
in the {Qun ^a) plane, which is necessary in order to try to 
reject the flatness hypothesis. This is shown below (Fig. 1, 
§3). 

Moreover, use of a linear regime comoving ruler also 
implies that peculiar velocity effects become negligible. 

Note that since the method uses data from a single sur- 
vey, it is qualitatively quite different from the concept of 
cosmic complementarity (Eisenstein, Hu & Tegmark 1998; 
Lineweaver 1998). 

The distance relations arc reviewed in §2, a method 
of illustrating the principle is explained in §3, results are 
presented in §4, and conclusions are made in §5. 

2. Distance relations 

Using the terminology of Weinberg (1972), the distance of 

use in the radial direction in comoving coordinates is the 
proper distance [eq. (14.2.21), Weinberg 1972; denoted Rx 
by Peebles (1993), eq. (13.28)], 



The curvature radius is then 



=-/' 

Hq Ji/ 



da 



(1) 



where c is the speed of light, iJo is the Hubble constant, 
Vlxa is the matter density divided by the critical density, 
Oa is the cosmological constant, z is redshift, and the 
dimensionless curvature is 



Rc 



1 



(3) 



For likely values of the curvature parameters, i.e. 
(Om ~ 0.3, Oa « 0.7), the value of Rc is constrained to 
Rc ^ 10h~^ Gpc for a curvature of <^ 0.1. Since 
the horizon radius, using proper distance d, is Rh ~ 
10h~^ Gpc, the present estimates of curvature can be suc- 
cintly restated as 



Rc^.Rh'^ lO/i-^ Gpc. 



(4) 



The proper distance is not usually useful in observa- 
tional cosmology, unless phenomena in comoving coordi- 
nates are being studied. The distances more commonly 

used are the proper motion distance [p. 485, Weinberg 
(1972); called 'angular size distance' by Peebles (1993), 
p.319, eq. (13.29)] 



Rcsinh[d{z)/Rc], Ok < 



dr,m{z) = { d{z), 



Ok = 



(5) 



Rcsm[d{z)/Rc], Ok>0. 



and distances which are greater or smaller than this by a 
factor of (1 -|- z). 

The tangential distance of use in comoving coordinate 
work can be written as 



d±{z,6e) = sedpuiiz) 

Rcsmh[d{z)/Rc], Ok < 
d{z), Ok = 

Rcsm[d{z)/Rc], Ok>0. 

where 66 is an angle in radians on the sky. 



(6) 



Ok = Om -|- Oa — 1- 



(2) 



3. Method 

Equations (1) and (6) clearly show that for a standard 
ruler placed at a large fraction of the curvature radius from 
the observer, the sin or sinh term will strongly distinguish 
the curved and flat cases. 

In order to test for independence of orientation of a 
standard ruler, consider a standard ruler of fixed comoving 
size L which is observationally detected near a redshift z, 
radially as a redshift interval 6z and tangentially as an 
angular size 66. Using equations (1) and (6), the length of 
the ruler placed in the radial direction is 

d{z + 6z,flm,i^A) - d{z,Slm,^A), (7) 
and the length in the tangential direction is 

d^{z, 66, nrr^QA), (8) 

where the implicit dependence on the curvature parame- 
ters is now made explicit. Define the difference between 

these as 

Ad(z,Om,OA)= (9) 
d±{z, 66, Om, Oa) - [d{z + 6z, Om, Oa) - d{z, Om, Oa)] . 
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Unless the Copernican principle is very surprisingly 
violated, Ad should be equal to zero, within the measure- 
ment uncertainties and intrinsic uncertainties of the ruler. 
If the Universe is 'slightly' curved, then a flat universe can 
be refuted if 







(10) 



can be refuted for all acceptable values of Om, or if the 
solutions to the equation are for values of d± inconsistent 
with L. 

The method is then defined by 

(a) choosing the acceptable range of flm, e.g. < ftm < 1, 

(b) choosing z, L and AL (the uncertainty in L), 

(c) calculating the implied values of 5z and 66, 

(d) assuming AL to be the total measurement and intrin- 
sic uncertainty (in length units) in each of the radial 
and tangential directions separately. 



and then calculating 

Ad{z,nm, 1 - f^m) 



(11) 



and the Gaussian probability that any of these values are 
consistent with zero, i.e. 'orientation- free'. 



Pr 



o-free 



erfc[|Ad|/(v^cr)], 



(12) 



where = 2(AL)^, for each value of fim- 

As is seen in Fig. 1 and discussed in §4, this is insuffi- 
cient on its own to rTilc out a flat universe hypothesis, so 

the additioual hypothesis that 

d±{z,5e,nm,l-^m)=L, (13) 

is required, where the Gaussian probability is 

it = evic[\d^{z, 66, Om, 1 - f^m) - L\/{V2a)], (14) 

whore cr^ = 2(AL)^ as above. (The radial ruler size could 
equally well be used here.) 

The independence of the ruler from orientation and the 
value of the length of the ruler are independent hypothe- 
ses, so the hypothesis of a flat universe is then rejected at 
the 

1 - ^^o-frcc^^L (15) 
confidence level. 

4. Results 

Fig. 1 shows that for reasonable values of (Om, ^^a) for 

a non-flat universe, the sin or sinh factor which relates 
the radial and tangential directions, making a ruler of dif- 
ferent sizes if the universe is curved [eq. 6], can be com- 
pensated for by the non-linearity in the distance redshift 
relation. In other words, even if the Universe really is 
'slightly' curved, a wrong pair of values (Omjf^A), where 
Qia is larger (smaller) than the true value for a hyperbolic 
(spherical) universe, can be found for which the radial and 




m 

Fig. 1. DiflFerences between tangential and radial lengths im- 
plied by assuming a flat universe, Ad{z, Clm, 1 — f2m) (eq. 10), 
and diflFerences between the tangential ruler size and its known 
size, d_L{z,5d,^lni, i — S^m) — L, both shown in comoving 
Mpc as a function of fim- The Ad curves increase with 
f2m; the d± ~ L curves decrease with fim- Various input hy- 
potheses for hyperbolic and spherical universe models and red- 
shifts z are shown with different curves as labelled. The points 
where the increasing curves pass through zero are those where 
a flat universe implies independence of the ruler size from ori- 
entation. The points where the decreasing curves pass through 
zero are those where the ruler size is the known (e.g. zero red- 
shift) size of the ruler. 



tangential sizes of the standard ruler are equal. These so- 
lutions are represented by the zero crossings of the curves 
which increase with f2ni in Fig. 1. The differences between 
the 'true' values and those required for independence 
of L from orientation (for a wrong, flat solution) are not 
large. 

This is why the size of the ruler needs to be known 
[eq. (13)]. The curves oi d± — L in Fig. 1 (which decrease 
with increasing Om) are reasonably steep near both the 
input and the 'orientation-free, flat' values of fini, but to 
distinguish the point of intersection of the curves from zero 
would require uncertainties of much less than 10 Mpc, 
or <C 10% of the ruler size. 

The probabilities for rejecting the flat universe based 
on the individual requirements of orientation indepen- 
dence [eq. (12)], correctness of the value of L [(eq. 14)] and 
the combined probabilities [eq. (15)] are shown in Figs 2 
and 3. In the latter figure, higher precision, i.e. 1%, is 
assumed in L. The minimum values of 1 — fg-free-'^L 
listed in Table 1. 

For an uncertainty of ~ 10%, none of the 'slightly' 
curved models, each having [$1^1 = 0.1 and Rc ~ Rh ~ 
10h~^ Gpc, would enable significant rejection of a fiat uni- 
verse. However, for an uncertainty of 1%, all three of the 
same 'slightly' curved models, using data at z = 3, would 
enable significant rejection of a flat universe, i.e. at the 
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Table 1. Minimum rejection values 1 — f'o-free-'^L "-"^ ^'-'^'^ curves in Figs 2 and 3. 



Hm 


nA 




z 

[h-^ Mpc) 


L 

(/i-i Mpc) 


AL 


min{l-Po.freeft} 


0.2 


0.7 


-0.1 


3 


130 


10 


31% 


0.4 


0.5 


-0.1 


3 


130 


10 


20% 


0.4 


0.7 


+0.1 


3 


130 


10 


24% 


0.4 


0.5 


-0.1 


2 


130 


10 


14% 


0.2 


0.7 


-0.1 


3 


130 


1.3 


99.7 % 


0.4 


0.5 


-0.1 


3 


130 


1.3 


94% 


0.4 


0.7 


+0.1 


3 


130 


1.3 


97% 


0.4 


0.5 


-0.1 


2 


130 


1.3 


80% 



Table 2. Minimum rejection values 1 — ^'o-free^' ^'^^ Table 1, but for a larger range in f2m, f^A, with z = 3, L = 130 h ^ Mpc, 
AL = 1.3 Mpc fixed. 



Om 




f2« 


min{l - Po-free^} 


fim 


Oa 


n« 


min{l - Po-free-PL) 








dependence 


on fim 






0.1 


0.8 


-0.1 


100.0% 


0.1 


1.0 


0.1 


100.0% 


0.2 


0.7 


-0.1 


99.7% 


0.2 


0.9 


0.1 


99.9% 


0.3 


0.6 


-0.1 


97.6% 


0.3 


0.8 


0.1 


99.2% 


0.4 


0.5 


-0.1 


93.9% 


0.4 


0.7 


0.1 


97.2% 


0.5 


0.4 


-0.1 


89.9% 


0.5 


0.6 


0.1 


93.8% 


0.6 


0.3 


-0.1 


86.4% 


0.6 


0.5 


0.1 


89.5% 


0.7 


0.2 


-0.1 


83.4% 


0.7 


0.4 


0.1 


84.9% 


0.8 


0.1 


-0.1 


77.5% 


0.8 


0.3 


0.1 


80.1% 


0.9 


0.0 


-0.1 


70.5% 


0.9 


0.2 


0.1 


75.3% 








a higher value of 








0.1 


0.7 


-0.2 


100.0% 


0.1 


1.1 


0.2 


100.0% 


0.2 


0.6 


-0.2 


100.0% 


0.2 


1.0 


0.2 


100.0% 


0.3 


0.5 


-0.2 


100.0% 


0.3 


0.9 


0.2 


100.0% 


0.4 


0.4 


-0.2 


99.9% 


0.4 


0.8 


0.2 


100.0% 


0.5 


0.3 


-0.2 


99.7% 


0.5 


0.7 


0.2 


100.0% 


0.6 


0.2 


-0.2 


99.2% 


0.6 


0.6 


0.2 


99.9% 


0.7 


0.1 


-0.2 


98.3% 


0.7 


0.5 


0.2 


99.6% 


0.8 


0.0 


-0.2 


97.2% 


0.8 


0.4 


0.2 


99.0% 


0.9 


-0.1 


-0.2 


95.6% 


0.9 


0.3 


0.2 


98.0% 



1 — P ^ 95% level. Lower redshift data {z = 2) only pro- 
vides marginally significant rejection. 

For the potentially most interesting case of AL = 1%L 
and z = 3, Table 2 shows how the strength of the rejection 
varies with and fi^- A nearly flat universe with a value 
of the matter density considerably higher than presently 
estimated, i.e. flra > 0.5, would be more difficult to re- 
ject, given a fixed absolute curvature Ifl^l = 0.1. On the 
other hand, a larger absolute curvature = 0.2) would 

enable refutation of a flat universe to better than 99.7% 
confidence for Om < 0.5. 



5. Discussion and Conclusions 

Does a standard ruler of the required precision exist? This 
depends, of course, on what the curvature of the Universe 
really is. If the Universe is 'exactly' flat in a geometrical 
sense, i.e. if the covering space is R^, then use of local 
geometrical techniques would not be sufficient to prove 
that the Universe is not curved. Proof that the Universe is 
flat and multiply connected (e.g. see Luminet & Roukema 



1999 for a review) would be one way of using global geom- 
etry to prove that the Universe is not (on average) curved. 

However, if the radius of curvature is no bigger than 
the horizon, then the calculations above show that for rea- 
sonable values of (f2m,fiA), a precision of 1% in the ap- 
plication of a standard ruler at a redshift of z = 3 would 
be sufficient. At higher redshifts, less precision would be 
needed, but the possibility of having large surveys includ- 
ing sufficient amounts of both radial and tangential stan- 
dard ruler information is unlikely in the next decade at 
2:> 3. 

At z = 3, surveys of quasars or of Lyman break galax- 
ies (e.g. Adelberger et al. 1998; Giavalisco et al. 1998) 
of sufficient quality to detect comoving features at large 
scales in the power spectrum should be feasible. 

Whether or not fixed, comoving features in the power 
spectrum of density perturbations, as traced by these ob- 
jects, exist and are detectable, is still a controversial sub- 
ject. Observations by several different groups suggest that 
a peak near the maximum in the power spectrum, at 
L « 130 ± 10 Mpc, is common to galaxies and super- 
clusters of galaxies at low redshift (Broadhurst et al. 1990; 
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Fig. 2. Rejection probabilities based on the results in Fig. 1, 
assuming that L = 130 Mpc, AL — 10 Mpc, and using 
equations (12), (14) and (15). Curve styles are as for Fig. 1. 
The thin curves are for 1 — Pq.^^^^ and 1 — P^, the bold curves 
are for 1 - Po-free-^L- 



I 




m 

Fig. 3. As for Fig. 2, for L = 130 h'^ Mpc, AL = 1.3 h'^ Mpc 

(i.e. 1% precision). The 1^^'o-frcc^ curves in bold are barely 
visible, close to the 100% rejection limit. Their minima are 
listed in Table 1. 



Broadhurst & JafFe 2000; da Costa ct al. 1993; Baugh 
& Efstathiou 1993; Gaztanaga & Baugh 1998; Einasto et 
al. 1994, 1997; Deng et al. 1996; Guzzo 1999; Tucker et 
al. 1998) and Lyman break objects (Broadhurst & JafFe 
2000) and quasars (Deng et al. 1994; Roukema & Mamon 
2000, 2001) at high redshift. Since 130 h'^ Mpc is above 
the present turnaround scale, it should be fixed in co- 
moving coordinates. Moreover, several possible theoreti- 
cal explanations for this feature, which would also imply 
other 'oscillations' in the power spectrum, include acous- 
tic oscillations in the baryon-photon fluid before last scat- 



Table 3. Statistical redshift accuracy Az needed for positions 
of low redshift structures used for calibrating L. For a given 
redshift z and given metric parameters, the uncertainty in red- 
shift Az corresponding to AL = 1.3 Mpc (1% precision in 
L = 130 h-^ Mpc) is given. This shows that over the domain 
shown, a precision of <,5xl0~* in redshift is required. 



^m ^^A ^ 

0.1 0.3 0.5 

0.3 0.0 0.0005 0.0006 0.0007 

0.3 0.7 0.0005 0.0005 0.0006 

1.0 0.0 0.0005 0.0006 0.0008 



tcring, in high baryon density models (Eiscnstcin et al. 
1998; Meiksin et al. 1998; Peebles 1999), and features from 
sub-Planck length physics which survive through to oscil- 
lations in the post-inflation power spectrum, for weakly 
coupled scalar field driven inflationary models (Martin & 
Brandenberger 2000a,b). 

No group yet claims that the precision of the peak 
is better than 10%. Observational improvements (homo- 
geneity of surveys, numbers of objects), refinements in 
statistical analysis techniques, and an unambiguous the- 
oretical explanation for the peak might help reduce the 
uncertainty in the value. 

What prospects for observational improvements from 
surveys expected to be completed within 1-5 years might 
potentially approach 1% precision in LI 

The best possible improvements in precision would 
presumably scale as Poisson errors, so that improving from 
~ 10% precision to 1% would require a factor of 100 in- 
crease in the numbers of objects relative to previous sur- 
veys. 

The continuation of the original observations by 
Broadhurst et al. (1990) which confirm the original re- 
sult include more than 1000 spectroscopic redshifts up to 
maximum redshifts of 2; ~ 0.4 (Broadhurst & JafFe 2000; 
Broadhurst 1999). The VIRMOS shallow survey (Le Fcvre 
et al. 2001) is expected to obtain 100,000 spectroscopic 
redshifts of galaxies in fields of size ~ 2° to a limiting 
magnitude of Iab ~ 22.5, i.e. with typical median (max- 
imum) redshifts of z ~ 0.6 {z ~ 1.3) (Crampton et al. 
1995). This is probably the best near future survey which 
more or less matches the conditions of the Broadhurst et 
al. 1990 observations. As long as the fact that the galax- 
ies will be spread over a somewhat larger redshift range 
does not adversely affect improvements in the estimation 
oF L, the factor of ~ 100 increase in numbers should be 
sufficient to provide the precision required. 

However, iF the cfFcct is anisotropic as suggested by 
Einasto et al. (1997), then the VIRMOS shallow survey 
may not be enough, since none of the four fields are close 
to the directions of the Broadhurst et al. flelds. 

Wide angle surveys may therefore be more useful. 

In the 2dF Galaxy Redshift Survey (2dFGRS, Colless 
et al. 1999) it is expected to observe spectroscopic red- 
shifts of 250,000 galaxies over 2000 sq.deg. with a mean 
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redshift oi z = 0.1. For the Sloan Digital Sky Survey 
(SDSS, Sloan et al. 2001), it is planned to observe 
about 1,000,000 galaxy spectroscopic rcdshifts over 10,000 
sq.deg. with a median redshift of z = 0.1. 

Both surveys clearly provide the increase in numbers 
of objects required. A possible problem in precision is 

the fact that only about 100 galaxics/sq.dcg. will be ob- 
served in these surveys as opposed to around 1000 galax- 
ies/sq.deg. in the Broadhurst et al. 1990 fields. If individ- 
ual structures ('walls', 'filaments') arc less sharply traced 
in position relative to the Broadhurst et al. 1990 surveys, 
then this could provide an additional noise factor. 

Will spectroscopic redshift accuracy be sufficient not to 
provide an additional source of uncertainty in the radial 
direction? Table 3 shows that a precision in redshift of 
Ax < 5x10""* corresponds to AL = 1.3 /i"^ Mpc. 

The VIRMOS survey (Le Fevre et al. 2001) is only 
expected to achieve spectral resolution of A/ A A ^ 250 
for the full 100,000 galaxies, implying Az - 4x10-^(1 + 
z)^10^^. Although a sub-sample will be observed at 
A/AA ^ 2500 — 5000 in order to study biases, it is hard 
to see how this can be sufficient for the purposes of ob- 
taining a precise value of L. It would be preferable if the 
full 100,000 galaxy sample could be observed at the higher 
resolution. 

On the other hand, the spectral resolutions of the 2dF- 
GRS and SDSS are expected to be A/AA ~ 1000 - 2000, 
so will provide approximately the precision required. 

However, Az<,5xl0~^ precision in redshift corre- 
sponds to Azc <, 300km/s. Typical galaxy velocity disper- 
sions in loose groups and clusters and typical bulk ve- 
locities, or in other words, typical galaxy velocities with 
respect to the comoving reference frame, are typically at 
about this scale or up to nearly an order of magnitude 
higher. So, averaging over these 'random' errors will be 
required in order to obtain the precision required. 

In the analysis of the low redshift surveys, it should be 
kept in mind that even though the redshifts are smaller 
than unity, a clear and precise detection of the L scale 
would require approximately correct values of the metric 
parameters, and may be missed if wrong values are used. 

It should also be noted that, in the hypothesis that the 
proposed ruler actually exists, and can be traced back to 
the primordial Universe, the theoretical explanation of the 
scale might heavily rely on some assumption on yet un- 
known or difficult to measure cosmological or fundamental 
physical parameters. This would introduce an additional 
uncertainty if the theoretical model were to be used to 
define the size of the standard ruler. 

In conclusion, if the Universe is indeed 'slightly' 
curved, then this method could potentially be promising 
for proving that the Universe is not flat. 

Note that small scale clustering on quasi-linear or non- 
linear scales (cf. Alcock & Paczynski 1979) would be more 
difficult to use as a standard ruler, due to evolution in the 
ruler length. 
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